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Abstract 

We discuss thermal and dynamical properties of Bose condensates confined by 
an external potential. First we analyze the Bose-Einstein transition tempera- 
ture for an ideal Bose gas in a generic power-law potential and d-dimensional 
space. Then we investigate the effect of the shape of the trapping potential 
on the properties of a weakly-interacting Bose condensate. We show that 
using exotic trapping potentials the condensate can exhibit interesting co- 
herent quantum phenomena, like superfluidity and tunneling. In particular, 
we consider toroidal and double-well potentials. The theoretical results are 
compared with recent experiments. 

1 Introduction 

In 1995 the Bose-Einstein condensation (BEC), i.e. the macroscopic occu- 
pation of the lowest single-particle state of a system of bosons, has been 
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experimentally achieved with clouds of confined alkali-metal atoms at ultra- 
low temperature (about 100 nK). In that year, three different groups (at 
JILA with 87 Rb atoms [1], at MIT with 23 Na atoms [2], and at Rice Univer- 
sity with 7 Li atoms [3]) have obtained the BEC by using the same technique: 
a laser cooling and confinement in a magnetic trap and an evaporative cool- 
ing. In 1997 the Nobel Prize in Physics has been given for the development 
of methods to cool and trap atoms with laser light [4]. Nowadays more than 
twenty experimental groups have achieved BEC by using different geometries 
of the confining trap and atomic species. 

The BEC phase transition has been predicted for a homogeneous ideal 
gas of Bosons by Einstein [5] on the basis of a paper of Bose [6]. The idea of 
Bose condensate has been used by London to describe the superfluid behavior 
of 4 He [7]. 4 He is a strongly-interacting system and the condensed fraction 
does not excedes 10%. Instead, in the dilute alkali-metal atoms of recent 
experiments, the condensed fraction can be more than 90%. Thus, dilute 
vapors are ideal systems to investigate thermal and dynamical properties of 
the Bose condensate and the role of the trapping potential. 

In this review paper we discuss some recent theoretical results we have 
obtained for Bose gases confined in potentials with exotic shapes, like power- 
law, toroidal and double-well potentials. Exotic potentials can be used in 
future experiments to find signatures of new coherent matter-wave phenom- 
ena. 
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2 Confined Bose gas and BEC temperature 



In this section we investigate a confined ideal quantum gas of Bosons. By 
using the semiclassical approximation, we derive analytical formulas for the 
BEC transition temperature. 

In the grand canonical ensemble of equilibrium statistical mechanics [8], 
the average number N a of Bosons in the single-particle state \a) with energy 
e a is given by 

1 

Na = e P(e a -n) _ i ' t 1 ) 
where \i is the chemical potential and (3 = l/(kT) with k the Boltzmann 
constant and T the absolute temperature. In general, given the single-particle 
function N a , the average total number N of particles of the system reads 

N = J2N a . (2) 

a 

This condition fixes the chemical potential /i. Thus, /j, is a function of (3 and 
N . fj, cannot be higher than the lowest single-particle energy level eo, i.e., it 
must be /i < e . When /i approaches e the function N becomes very large 
and consequently iV becomes of the same order of N. The physical meaning 
is that the lowest single-particle state becomes macroscopically occupied and 
one has the so-called Bose-Einstein condensation (BEC). It is clear that if 
\x = e then No diverges. It is a standard procedure to calculate the condensed 
fraction N /N and also the BEC transition temperature Tb by studying the 
non divergent quantity iV — N at ji = eo as a function of the temperature. 
This procedure is particularly effective by using the semiclassical approxima- 
tion. 
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In the semiclassical approximation, the system is described by a con- 
tinuum of states and, instead of e a , one uses the classical single-particle 
phase-space energy 

^,p) = |Ut/(r), (3) 

where p 2 /(2m) is the kinetic energy and U(r) is the confining external po- 
tential. In this way one obtains the semiclassical single-particle phase-space 
distribution of Bosons 

w ( r » P) = e/^.p)-,) _ x • ( 4 ) 

Note that the accuracy of the semiclassical approximation is expected to be 
good if the number of particles is large and the energy level spacing is smaller 
then kT [9]. 

For the sake of generality, let us consider a d-dimensional space. Then 
the quantum elementary volume of the single-particle 2d-dimensional phase- 
space is given by (2Ti%) d ) where U is the Planck constant [9]. It follows that 
the average number N of Bosons in the d-dimensional space can be written 

as 



N 
where 



d d p 



is the spatial distribution, A = (2irh P/m) 1 / 2 is the thermal length, 
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r(n) 

is the Bose function, and is the factorial function. 



It is important to observe that, while Eq. (4) is divergent at /i = e(r, p), 
Eq. (6) is not divergent at /j, — U(r) because gd/2^) = ((d/2), where ((x) 
is the Riemann zeta-function. It means that the integration over momenta 
removes the divergent contribtion coming from the ground-state. It follows 
that the Eq. (5) describes only the total number of non-condensed Bosons. 
Note that this number can also be written as 



where p(e) is the density of states. It can be obtained from the semiclassical 
formula 




where 5(x) is the Dirac delta function and T(n) is the factorial function. 
This result is the generalization of the formula for ideal homogeneous Bose 
gases in a box of volume V [8]. It shows that, in the semiclassical limit, the 
non- homogeneous formula is obtained with the substitution ji — > /i — U(r), 
also called local density approximation. 

2.1 BEC temperature for power- law potentials 

In many experiments with alkali-metal atoms, the external trap can be accu- 
rately described by a harmonic potential. More generally, one can consider 
power-law potentials, given by 

U(r) = Ar n = -(^ 
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where n is the power-law exponent and A is the trap constant. Clearly, 
with n = 2 one gets the harmonic potential. Here we have introduced an 
energy parameter hcu and a length parameter r , which can be chosen as 



r = \Jh/ (muo). The power-law potential is interesting for studying the 
effects of adiabatic changes in the trap. The density of states of a quantum 
gas in the power-law potential can be calculated from Eq. (8) and reads 

p(e) = - d '^^P „ (^o)^ • (10) 

As previously stated, Eq. (5) and (7) give the total number of Bosons 
only above critical temperature T B . Below T B , Eq. (5) and (7) describe 
the non-condensed number of particles N — N , where N is the number of 
condensed particles. At the Bose transition temperature T B , the chemical 
potential is /i — and in this way one obtains the critical temperature. In 
the case of the power-law potential the critical temperature is 

2n 

kT B = c{d,n) huo N ' d (™+ 2 ) , (11) 
where c(d, n) is a numerical coefficient given by 

2#r(f + 1) 



c(d, n) 



d(n + 2) 



(12) 



Below the critical temperature T B , one has N ^ and from Eq. (5) and (7) 
one gets the T dependence of the condensed fraction 

d(n+2) 

N \T B ) 
where N is the number of Bosons in the gas. 



It is important to observe that from the previous formulas one easily 
derives the thermodynamic properties of quantum gases in harmonic traps 
and in a rigid box. In fact, by setting n = 2 one gets the formulas for the Bose 
and Fermi gases in a harmonic trap (in the case of a anisotropic harmonic 
potential, ujq is the geometric average of the frequencies of the trap). The 
results for a rigid box are instead obtained by letting - — > 0, where the 
density of particles per unit length is given by N/VL d and Vt d = /T(| + 1) 
is the volume of the d-dimensional unit sphere. 

Finally, one notes that in the formula of the BEC transition temperature 
Tb it appears the function £(| + -). Because < oo for x > 1 but 

C(l) = oo, one easily deduces the following theorem [9]. 

Theorem. Let us consider an ideal Bose gas in a power-law isotropic 
potential U(r) = A r n with r = |r| = (J2f=i x \ ) 1 ^ 2 - BEC is possible if and 
only if the following condition is satisfied 

d{n + 2) 

where d is the space dimension and n is the exponent of the confining power- 
law potential. 

This is a remarkable inequality. For example, for d — 2 one finds the 
familiar result that there is no BEC in a homogeneous gas — > 0) but BEC 
is possible in a harmonic trap (n — 2). Moreover, one obtains that for d — 1 
BEC is possible with 1 < n < 2. 
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3 Bose condensate in external potential 

In this section we study in detail the static and dynamical properties of a 
trapped dilute Bose condensate at zero temperature. Note that, if the system 
of Bosons is dilute, at zero temperature the non-condensed fraction can be 
neglected [10]. 

Let us consider a iV-body quantum system with Hamiltonian H . The 
exact time-dependent Schrodinger equation can be obtained by imposing the 
quantum least action principle to the action 

r d 

S = J dt< ^{t)\ih— - H\$(t) > , (13) 

where $ is the many-body wave-function of the system and 

< >= J d 3 r 1 ...d 3 r N $*(r 1 ,...,r N ,t)A<S>(r 1 ,...,r N ,t) , (14) 

for any quantum operator A. Looking for stationary points of S with respect 
to variation of the conjugate wave-function <3>* gives 

d 

ih— $(n,...,rjv,f) = H $(ri,...,rjv,f) , (15) 

which is the many-body time-dependent Schrodinger equation. 

As is well known, except for integrable systems, it is impossible to obtain 
the exact solution of the many-body Schrodinger equation and some ap- 
proximation must be used. Here we discuss the zero-temperature mean-field 
approximation for a system of trapped weakly-interacting bosons in the same 
quantum state, i.e. a Bose-Einstein condensate. In this case the Hartree-Fock 
equations reduce to only one equation, the Gross-Pitaevskii equation, which 
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describes the dynamics of the condensate. As previously discussed, this equa- 
tion is intensively studied because of the recent experimental achievement of 
Bose-Einstein condensation for atomic gasses in magnetic traps at very low 
temperatures. 

The Hamiltonian operator of a system of N identical bosons of mass m 
is given by 

N t2 i N 

H = E ( - ^\ + U(Ti)) + - E V(r it r,) , (16) 

i=l ij=l 

where U(r) is an external potential and V(r, r') is the interaction potential. 
In the mean-field approximation the totally symmetric many-particle wave- 
function of the Bose-Einstein condensate reads 

$(ri, .., r N , t) = *(n, t)*(r 2 , t) ... *( rjv _i, t)V(r N , t) , (17) 

where ^(r,t) is the single-particle wave-function. Note that such factor- 
ization of the total wave-function is exact in the case of a non-interacting 
condensate. The quantum action of the system is then simply given by 

S GP = N Jdt< *{t)\ih— - h a \*(t) > , (18) 

where 

h s = -^_V 2 + U{r)+ l -{N-l) jd 3 r'\^(v',t)\ 2 V(r,r'). (19) 

We call Sgp the Gross-Pitaevskii (GP) action of the Bose condensate [10]. By 
using the quantum least action principle we get the Euler-Lagrange equation 

ih^-V(r,t) = [- |-V 2 + U(r) + (N - 1) f dW(r, r')|*(r', t)| 2 ]tf(r, t) , 

(20) 



which is an integro-differential nonlinear Schrodinger equation. Such equa- 
tion and the effect of a finite-range interaction have been analyzed only by 
few authors [11,12]. In fact, at low energies, it is possible to substitute the 
true interaction with a pseudo-potential 

V(r,r>)=g5 3 (r-r>), (21) 

where g = 4irh 2 a s /m is the scattering amplitude and a s the s-wave scattering 
length. The scattering length is positive (repulsive interaction) for 87 Rb and 
23 Na atoms but negative (attractive interaction) for 7 Li atoms. Moreover, 
for large N, the factor (N — 1) can be substituted with N. In this way one 
obtains the so-called time-dependent GP equation 



d 
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iKft*(r, t)=[- — V 2 + U(r) + gN\*(r, t)\ 2 \ tt(r, t) , (22) 

that is the starting point of many calculations [13]. Note that the GP equa- 
tion is accurate to describe the condensate of weakly-interacting bosons only 
near zero temperature, where thermal excitations can be neglected. The 
ground-state solution ip(r) of the GP equation is found setting \I/(r, t) = 
e -W^( r ) j n the previous equation. In this way one gets the stationary GP 
equation 

■ " 2^ V2 + C/(r) + 9N\^)\ 2 ]^) = V V>(r) , (23) 
where \i is the chemical potential. Note that this equation can be also found 
by mimimizing the GP energy functional 

E = N Jd 3 r |^|VV(r)| 2 + U(rM(v)\ 2 + ^iV|^(r)| 4 , (24) 
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with the normalization condition 

|rf 3 r|^(r)| 2 = l, (25) 

which fixes the chemical potential /i, that is the Lagrange multiplier of the 
stationary minimization problem. 

To calculate the energy and wavefunction of the zero-temperature elemen- 
tary excitations, one must solve the so-called Bogoliubov-de Gennes (BdG) 
equations [14]. The BdG equations can be obtained from the linearized time- 
dependent GP equation. Namely, one can look for zero angular momentum 
solutions of the form 

#(r, t) = e~>* [V>(r) + u(r)e _ihrt + v*{v)e itjjt ] , (26) 

corresponding to small oscillations of the wavefunction around the ground 
state solution ip(r). By keeping terms linear in the complex functions u(r) 
and v(r), one finds the following BdG equations 

r h 2 

- — V 2 + U(p, z)-n + 2gNWr)\ 2 u(r) + gN\i>(r)\ 2 v(r) = hw u(r) , 
L 2m J 

r h 2 

V 2 + f/(p,^)- / u+2c/iV|^(r)| 2 w(r)+^iV|^(r)| 2 M(r) = -hw v(r) . (27) 

L 2m J 

The BdG equations allow one to calculate the eigenfrequencies ui and hence 
the energies Hou of the elementary excitations. This procedure is equiva- 
lent to the diagonalization of the N-body Hamiltonian of the system in the 
Bogoliubov approximation. 
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3.1 Bose condensate in a toroidal potential 

In this section we consider a 3-D toroidal trap given by a quartic Mexican 
hat potential along the cylindrical radius and a harmonic potential along 
the z axis. As we have recently shown [15], the resulting toroidal trapping 
potential is very flexible and it is possible to modify considerably the density 
profile of the BEC by changing the parameters of the potential or the number 
of bosons. Moreover, the toroidal trap can be used to create a superfluid, 
namely persistent currents in absence of imposed rotation. 

The toroidal trap we are discussing can be described by the following 
potential 

U(v) = iA(p 2 - pi? + imc^V , (28) 



where p = \/x 2 + y 2 and z are the cylindrical coordinates. The potential 
U (r) is minimum along the circle of radius p = p at z = and it has a local 
maximum at the origin in the (x, y) plane. Small oscillations in the (x, y) 
plane around p have a frequency uj±_ = p (2\/m) 1 ^ 2 . 

In the Thomas-Fermi approximation, i.e. neglecting the kinetic energy in 
the stationary GP equation, we find 



1 



V2, 



^ (r) = ^iv (/U ~ ^ (r)) J 0(/U ~ ^ (r)) ' (29) 

where Q(x) is the step function. For our system we obtain that: a) the wave 
function has its maximum value at p = p and z = 0; b) for p < Apg/4 the 
wave function has a toroidal shape; c) for p > Apg/4 the wave function has 
a local minimum at p = z = 0; d) the chemical potential scales as p ~ N 1 ^ 2 . 
It is important to note that the TF approximation neglects tunneling effects: 
to include these processes, it is necessary to analyze the full GP problem. 
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We have performed the numerical minimization of the GP functional by 
using the steepest descent method. This method consists of projecting onto 
the minimum of the energy functional an initial trial state by propagating 
it in imaginary time. We have discretize the space with a grid of points 
taking advantage of the cylindrical symmetry of the problem. We have used 
grids up to 200 x 200 points verifying that the results do not depend on 
the discretization parameters. The number of iterations in imaginary time 
depends on the degree of convergence required and the goodness of the initial 
trial wave function. In calculations we have adopted ^-harmonic oscillator 
units. We write po m units a z = (/^/(mcc^)) 1 / 2 = 1 /im, A in units (hco z )a~ 4 = 
0.477 (5.92) peV/pm 4 and the energy in units hu z = 0.477 (5.92) peV for 
87 Rb ( 7 Li). Moreover, we have used the following values for the scattering 
length: a s = 50 (-13) A for 87 Rb ( 7 Li). 

In the case of positive scattering length ( 87 Rb) we can control the density 
profile of the BEC by modifying the parameters of the potential and also the 
number of particles. For small number of particles the condensate is essen- 
tially confined along the minimum of U(r), there is a very small probability 
of finding particles in the center of the trap so that the system is effectively 
multiply connected. As N increases the center of the trap starts to fill up 
and the system becomes simply connected. The value of N for which there 
is a crossover between the two regimes increases with the value of A and of 
Po and, within Thomas- Fermi approximation, scales like A 3//2 Pq. 

In the case of negative scattering length ( 7 Li), it is well known that for the 
BEC in harmonic potential there is a critical number of bosons N c , beyond 
which there is the collapse of the wave function (see also [16]). We obtain 

13 



the same qualitative behavior for the 7 Li condensate in our Mexican hat 
potential. However, in cylindrical symmetry, the collapse occurs along the 
line which characterizes the minima of the external potential, i.e. at p = p 
and z — 0. We notice that, for a fixed p , the critical number of bosons N c 
is only weakly dependent on the height of the barrier of the Mexican poten- 
tial. These results suggest that we can not use toroidal traps to significantly 
enhance the metastability of the BEC with negative scattering length. 

It is useful to study states having a vortex line along the z axis and all 
bosons flowing around it with quantized circulation. The observation of these 
vortex states is a signature of macroscopic phase coherence of trapped BEC. 
The axially symmetric condensate wave function can be written as 

M*)=MP,z)e iM , (30) 

where 9 is the angle around the z axis and k is the integer quantum number 
of circulation. The resulting GP energy functional can be written in terms 
of ^fc(r) by taking advantage of the cylindrical symmetry of the problem: 



+ 



h 2 T,dif>k(p,z) |2 dip k {p,z) 2 
F 2m L 1 

h 2 k 2 



E = N J pdpdzdO 



dp 1 q z 



+ 



Imp 



2 



+ U(p,z) \Mp,z)\ 2 + t;9N\Mp,z)\ 4 ■ (31) 



2 



Due to the presence of the centrifugal term, the solution of this equation for 
k 7^ has to vanish on the z axis providing a signature of the vortex state. 

As previously stated, vortex states are important to characterize the 
macroscopic quantum phase coherence and also superfluid properties of Bose 
systems. It is not difficult to calculate the critical frequency Q c at which a 
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vortex can be produced. One has to compare the energy of a vortex state 
in a frame rotating with angular frequency Q, that is E — QL Z , with the 
energy of the ground state with no vortices. Since the angular momentum 
per particle is hk, the critical frequency is given by Ml c = (E k /N — E /N)/k, 
where E^/N is the energy per particle of the vortex with quantum number k. 
In Table 1 we show some numerical results for vortices of 87 Rb. The critical 
frequency turns out to increase slightly with the number of atoms. This cor- 
responds to a moderate lowering of the momentum of inertia per unit mass 
of the condensate when N grows. 



TV 


EJN 




hfi c 


5000 


6.00 


7.87 


0.15 


10000 


7.61 


10.44 


0.16 


20000 


10.02 


14.22 


0.18 


30000 


11.93 


17.20 


0.20 


40000 


13.57 


19.75 


0.21 


50000 


15.04 


22.04 


0.23 



Table 1. Toroidal trap. Vortex states of 87 Rb atoms with k — 1 in the 
toroidal trap with p = 2 and A = 4 within Hartree-Fock approximation. 
Chemical potential and energy are in units of %oo z = 0.477 peV (oo z = 0.729 
kHz). Lengths are in units of a z — 1 /im. 

Once a vortex has been produced, the Bose condensate is superfluid if 
the circulating flow persists, in a metastable state, in the absence of an 
externally imposed rotation. Vortex solutions centered in harmonic traps 
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have been found, but such states turn out to be unstable to single particle 
excitations out of the condensate. The vortex state is superfluid if only 
if its lowest elementary excitation is positive [17]. To have the complete 
spectrum of elementary excitations, one must solve the BdG equations (see 
Eq. (27)). We have solved the two BdG eigenvalue equations by finite- 
difference discretization with a lattice of 40 x 40 points in the (p, z) plane. In 
this way, the eigenvalue problem reduces to the diagonalization of a 3200 x 
3200 real matrix. 



N 


5 x 10 3 


10 4 


2 x 10 4 


3 x 10 4 4 x 10 4 


5 x 10 4 


Tiuj 


1.22 


1.48 


1.73 


1.88 1.99 


2.08 



Table 2. Toroidal trap. Bogoliubov elementary excitation for a vortex state 
of 87 Rb atoms with k — 1 in the toroidal trap with p = 2 and A = 4. Units 
as in Tab. 1. 

We have tested our program in simple models by comparing numerical 
results with the analytical solution and verified that a 40 x 40 mesh already 
gives reliable results for the lowest part of the spectrum. The results are 
shown in Table 2: The lowest Bogoliubov excitation is always positive. We 
have also verified that vortex states become unstable by strongly reducing 
either density (down to about one hundred bosons in our model trap) or 
scattering length. Thus, we can conclude that a strongly-interacting Bose 
condensate in a toroidal potential is superfluid. 
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3.2 Bose condensate in a double-well potential 

In a recent experiment at MIT [18], the macroscopic interference of two Bose 
condensates released from the double minimum potential has been demon- 
strated. Here, we consider the same double-well potential. We analyze the 
ground state properties of the condensate and calculate the spectrum of the 
Bogoliubov elementary excitations as a function of the double-well barrier. 
By varying the strength of the barrier one can observe macroscopic quantum 
effects, like the formation of two Bose condensates, the collective oscillations 
and the quantum tunneling [19]. 

The double-well trap is given by a harmonic anisotropic potential plus a 
Gaussian barrier along the z axis, which models the effect of a laser beam 
perpendicular to the long axis of the condensate: 



where p = \Jx 2 + y' 2 , z and the angle 8 are the cylindrical coordinates. The 
parameter values appropriate for Ref. [18] are u p = 2tt x 250 Hz, oj z = 2n x 19 
Hz, and a = 6 fim. The anisotropic harmonic trap implies a cigar-shaped 
condensate (A = uo z /oo p = 15/250 < 1), where z is the long axis, and the 
Gaussian barrier of strength Uq creates a double-well potential. 

We have performed the numerical minimization of the GP energy func- 
tional by using the steepest descent method. At each time step the matrix 
elements entering the Hamiltonian are evaluated by means of finite-difference 
approximants using a grid of 200 x 800 points. In our calculations we have 
used the z-harmonic oscillator units. For 23 Na atoms, the harmonic length is 
a z = (^/(mwz)) 1 ' 2 = 4.63 /iin and the energy is hu z = 0.78 peV. Moreover, 




(32) 
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we have used the following value for the scattering length: a s = 3 nm [3]. 
Most of our computations has been performed for N = 5 x 10 6 atoms, a value 
typical of the MIT experiment [19]. 

In Figure 1 we show the ground state density profile of the 23 Na con- 
densate for different values of the strength of the barrier. By increasing the 
strength, the fraction of 23 Na atoms decreases in the central region and the 
Bose condensate separates in two condensates. The numerically calculated 
density profiles are in good agreement with the phase-contrast images of the 
MIT experiment [18] and with the Thomas- Fermi (TF) approximation, which 
neglects the kinetic term in the GP equation. Due to the large number of 
atoms involved (N = 5 x 10 6 ), only near the borders of the wave function 
there are small deviations from the TF approximation. Note that the poten- 
tial barrier U can be written as U c /k B = (37 /j,K)P/ a 2 (/zm 2 /mW), where 
P is the total power of the laser beam perpendicular to the long axis of the 
condensate and o = 6/zm is the beam radius. The conversion factor is 0.09 
mW/ (hu z ), such that U = 100 (in Hcu z units) gives a laser Power P = 9 mW 
[18]. 

Another important property of the BEC is the spectrum of elementary 
excitations. To calculate the energy and wavefunction of elementary excita- 
tions, one must solve the BdG equations. The excitations can be classified 
according to their parity with respect to the symmetry z — > —z. We have 
solved the two BdG eigenvalue equations by finite-difference discretization 
with a lattice of 40 x 40 points in the (p, z) plane. In this way, the eigenvalue 
problem reduces to the diagonalization of a 3200 x 3200 real matrix. We 
have tested our program in simple models by comparing numerical results 
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Figure 1: Double- well trap. Particle probability density in the ground state 
of N = 5 x 10 6 23 Na atoms as a function of the z axis at r = (symmetry 
plane). The curves correspond to increasing values of the strength U of 
the barrier (from to 500), in units of hu z = 0.78 peV. The laser power is 
given by the conversion formula P = 0.09 x Uq mW. Lengths are in units of 
a z = 4.63 fim. 
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Figure 2: Double-well trap. Lowest elementary excitation hui vs barrier 
energy U for N = 5 x 10 6 23 Na atoms. Energies in units of %oo z = 0.78 peV. 
The laser power is given by the conversion formula P = 0.09 x U mW. 

with analytical solutions [19]. 

In Figure 2 we show the lowest elementary excitations of the Bogoliubov 
spectrum for the ground state of the system. When the Gaussian barrier 
is switched off, one observes the presence of an odd excitation at energy 
quite close to Tiu = 1 (in units hu z ). This mode is related to the oscillation 
of the center of mass of the condensate along the z-axis, due to the har- 
monic confinement. This collective oscillation is an exact eigenmode of the 
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problem characterized by the frequency u z , independently of the strength 
of the interaction. The inclusion of the Gaussian barrier modifies the har- 
monic confinement along the z-axis and this odd collective mode decreases 
by increasing the strength of the barrier. 

For large values of the Gaussian barrier, i.e. when the BEC separates 
in two condensates, we find quasi-degenerate pairs of elementary excitations 
(even-odd). The lowest elementary excitation and the ground-state of the 
GP equation constitute one of such pairs and get closer and closer as the 
barrier is increased. This is not surprising because in the infinite barrier 
limit we have two equal and independent Bose condensates with the same 
energy spectrum. 

An interesting aspect of BEC in double-well traps is the possibility to de- 
tect the macroscopic quantum tunneling (MQT). The MQT has been recently 
investigated by Smerzi et al. [20]. They have found that the time- dependent 
behavior of the condensate in the tunneling energy range can be described 
by the two-mode equations 

z = —VI — z 2 sin0 , = Az H — cos , (33) 

V 1 — z 2 

where z = (Aq — N 2 )/N is the fractional population imbalance of the con- 
densate in the two wells, = 0i — 02 is the relative phase (which can be 
initially zero), and A = 4E mt / AE° . E mt is the interaction energy of the 
condensate and AE° is the kinetic+potential energy splitting between the 
ground state and the quasi-degenerate odd first excited state of the GP equa- 
tion. For a fixed A (A > 2), there exists a critical z c = 2V 'A — 1/A such that 
for < z « z c there are Josephson-like oscillations of the condensate with 

21 



period r = To/y/l + A, where r = 2nh/AE°. But for z c < z < 1 there is 
macroscopic quantum self-trapping (MQST) of the condensate: even if the 
populations of the two wells are initially set in an asymmetric state (z 7^ 0) 
they maintain the original population imbalance without transferring parti- 
cles through the barrier as expected for a free Bose gas. 



a s /af a 


A 


To (sec) 


Z c 


10- 1 


1108.337 


14.583 


0.060 


io- 2 


133.643 


13.887 


0.173 


10- 3 


1.390 


13.842 


none 


10" 4 


0.103 


10.253 


none 



Table 3. Double-well trap. Parameters of the MQT for different values 
of the scattering length a s with a^ a = 3 nm and t = 2iTh/AE . Condensate 
with N = 5 x 10 3 atoms. Barrier with Uq = 20 and a = 1.5 /jm. Energy 
barrier U in units of hu z = 0.78 peV. 

By solving the stationary GP equation in the MIT double-well trap with 
23 Na, we have found that the parameter A is larger than 10 4 also when few 
particles are present. Nevertheless, one can control the dynamics of the con- 
densate by reducing the scattering length a s and the thickness a of the laser 
beam. In particular, as shown in Table 3, the parameter A scales linearly 
with a s . This is an important point because recently it was confirmed exper- 
imentally the fact that it is now possible to control the two-body scattering 
length by placing atoms in an external field. This fact opens the way to a 
direct observation of a macroscopic quantum tunneling of thousands of atoms 
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through a potential barrier. 

4 Conclusions 

We have studied Bose gases in various trapping potentials above and below 
the BEC transition temperature. By using the semiclassical approximation, 
we have derived analytical formulas for the BEC transition temperature of 
an ideal Bose gas in a generic power-law potential and d-dimensional space. 
Then we have analyzed the effects of the shape of the external potential on 
a zero-temperature Bose condensate by using the Gross-Pitaevskii equation, 
which describes the macroscopic wave-function of the condensate, and the 
Bogoliubov-de Gennnes equations, which describe the elementary excitations 
of the condensate. We have shown that, contrary to the case of a Bose 
condensate in a harmonic potential, a strongly-interacting condensate in a 
toroidal potential is superfluid. Finally, we have investigated the conditions 
under which it is possible to find macroscopic quantum tunneling with a Bose 
condensate in a double-well potential. 

References 

[1] M.H. Anderson, J.R. Ensher, M.R. Matthews, C.E. Wieman, and E.A. 
Cornell, Science 269, 189 (1995). 

[2] K.B. Davis, M.O. Mewes, M.R. Andrews, N.J. van Druten, D.S. Drufee, 
D.M. Kurn, and W. Ketterle, Phys. Rev. Lett. 75, 3969 (1995). 



23 



[3] C.C. Bradley, C.A. Sackett, J.J. Toilet, and R.G. Hulet, Phys. Rev. 
Lett. 75, 1687 (1995). 

[4] S. Chu, C. Cohen-Tannouji, and W.D. Phillips, Nobel Prize in Physics 

(1997) . 

[5] A. Einstein, Preussische Akademie der Wissenshaften 22 261 (1924); 1 
3 (1925); 3 18 (1925) 

[6] S.N. Bose, Z. Phys. 26 178 (1924) 

[7] F. London, Nature 141 643 (1938); F. London, Phys. Rev. 54 947 
(1938); F. London J. Phys. Chem. 43 49 (1938). 

[8] K. Huang, Statistical Mechanics (John Wiley, New York, 1987); A.L. 
Fetter and J.D. Walecka, Quantum Theory of Many-Particle Systems 
(Mc Graw-Hill, Boston, 1971) 

[9] L. Salasnich, J. Math. Phys. 41 8016 (2000). 

[10] L. Salasnich, Int. J. Mod. Phys. B 14 1 (2000). 

[11] A. Parola, L. Salasnich and L. Reatto, Phys. Rev. A 57, R3180 (1998); 
L. Reatto, A. Parola and L. Salasnich, J. Low Temp. Phys. 113, 195 

(1998) . 

[12] L. Salasnich, Phys. Rev. A 61, 015601 (2000). 

[13] E.P Gross, Nuovo Cimento 20, 454 (1961); L.P Pitaevskii, Sov. Phys. 
JETP 13, 451 (1961). 



24 



[14] N.N. Bogoliubov, J. Phys. U.S.S.R. 11 23 (1947). 

[15] L. Salasnich, A. Parola and L. Reatto, Phys. Rev. A 59, 2990 (1999). 

[16] L. Salasnich, Mod. Phys. Lett. B 11 1249 (1997); L. Salasnich, Mod. 
Phys. Lett. B 12, 649 (1998). 

[17] D.S. Rokhsar, Phys. Rev. Lett. 79, 2164 (1997). 

[18] M.R. Andrews, C.G. Townsend, H.J. Miesner, D.S. Drufee, D.M. Kurn, 
and W. Ketterle, Science 275, 637 (1997). 

[19] L. Salasnich, A. Parola and L. Reatto, Phys. Rev. A 60 4171 (1999). 

[20] A. Smerzi, S. Fantoni, S. Giovannazzi and S.R. Shenoy, Phys. Rev. 
Lett. 79, 4950 (1997). 



25 



